Parallel Solution of Narrow Banded Diagonally Dominant Linear Systems
نویسندگان
چکیده
Abstract ScaLAPACK contains a pair of routines for solving systems which are narrow banded and diagonally dominant by rows. Mathematically, the algorithm is block cyclic reduction. The ScaLAPACK implementation can be improved using incomplete, rather than complete block cyclic reduction. If the matrix is strictly dominant by rows, then the truncation error can be bounded directly in terms of the dominance factor and the size of the partitions. Our analysis includes new results applicable in our ongoing work of developing an efficient parallel solver.
منابع مشابه
The Stable Parallel Solution of Narrow Banded Linear Systems
We report on experiments with a stable algorithm for the parallel solution of banded and periodically banded linear systems. In contrast to most of the known parallel algorithms which are stable only for symmetric positive deenite or diagonally dominant systems, the new algorithm incorporates pivoting without sacriicing eeciency. We report on numerical experiments conducted on an Intel Paragon,...
متن کاملA Comparison of Parallel Solvers for Diagonally Dominant and General Narrow-Banded Linear Systems II
We continue the comparison of parallel algorithms for solving diagonally dominant and general narrow-banded linear systems of equations that we started in [2]. The solvers compared are the banded system solvers of ScaLAPACK [6] and those investigated by Arbenz and Hegland [1, 5]. We present the numerical experiments that we conducted on the IBM SP/2.
متن کاملA Comparison of Parallel Solvers for Diagonally Dominant and General Narrow-Banded Linear Systems
We investigate and compare stable parallel algorithms for solving diagonally dominant and general narrow-banded linear systems of equations. Narrow-banded means that the bandwidth is very small compared with the matrix order and is typically between 1 and 100. The solvers compared are the banded system solvers of ScaLAPACK [11] and those investigated by Arbenz and Hegland [3, 6]. For the diagon...
متن کاملThe Explicit Spike Algorithm: Iterative Solution of the Reduced System
The explicit SPIKE algorithm applies to narrow banded linear systems which are strictly diagonally dominant by rows. The parallel bottleneck is the solution of the so-called reduced system which is block tridiagonal and strictly diagonally dominant by rows. The reduced system can be solved iteratively using the truncated reduced system matrix as a preconditioner. In this paper we derive a tight...
متن کاملThe stable parallel solution of general narrow banded linear systems
A list of technical reports, including some abstracts and copies of some full reports may be found at: A parallel architecture for query processing over a terabyte of text. June 1996. Abstract. We propose a stable algorithm for the parallel solution of banded and periodically banded linear systems. While most of the known parallel algorithms are stable only for symmetric positive deenite or dia...
متن کامل